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ABSTRACT 

Nitric acid is a possible biomarker in the atmospheres of exoplanets. An accurate 
line list of rotational and rotational-vibrational transitions is computed for nitric acid 
(HNO 3 ). This line list covers wavelengths longer than 1.42 /rm (0 - 7000 cm“^) and 
temperatures up to 500 K. The line list is computed using a hybrid variational - 
perturbation theory and empirically tuned potential energy and dipole surfaces. It 
comprises almost 7 billion transitions involving rotations up to J = 100. Comparisons 
with spectra from the HITRAN and PNNL databases demonstrate the accuracy of 
our calculations. Synthetic spectra of water - nitric acid mixtures suggest that nitric 
acid has features at 7.5 and 11.25 /rm that are capable of providing a clear signature 
for HNO 3 ; the feature at 11.25 /rm is particularly promising. Partition functions plus 
full line lists of transitions are made available in an electronic form as supplementary 
data to the article and at www. exomol. com. 

Key words: molecular data; opacity; astronomical data bases: miscellaneous; planets 
and satellites: atmospheres. 


1 INTRODUCTION 

The infrared spectrum of the nitric acid molecule, HNO 3 , 
is of astrophysical interest because of the growing interest 
in the study the of atmospheres of extrasolar planets. In 
particular, this interest is particularly linked to that of the 
planets with earth-like atmospheres with a high content of 
nitrogen and oxygen. Nitric acid ca n be clearly observed in 
the earth’s atmosphere from space (ICooper et ahlkoill l as 
it has a number of features which lie in water transparency 
windows. Its spectrum constitutes a possible biomarker as 
its existence is indicative of the presence of free oxygen and 
nitrogen, both of which lack strong spectral signatures at 
long wavelengths. It is also thought likely to pay an impor¬ 
tant rol e in nitrogen fixation on Mars and other earth-like 
planets JSummers fc Kharel2007| l . It has also been suggested 
that 7 -ray bursts could produce large q uantities of HNO 3 
in th e atmospheres of earth-like planets (iThomas fc Melot^ 
kOOd l. Nitric acid is a constituent of atmospheric ices on 
Earth and is thought likely t o be present in the ice crust 
of Europa (ICooper et al.ll^llll . Its for mation in interstella r 
ammonia ices has also been suggested JZanchet et aLllioiJ) , 
although nitric acid has yet to be observed in the interstellar 
medium. 

Given the importance of HNO 3 in the earth’s atmo¬ 
sphere, its infrared spectrum has been well-studied in the 
laboratory. However, many of these spectra remain incom¬ 
pletely analysed with very few assignments to transitions 
at wavelengths shorter than 5 pm. This lack of assignment 


means that the spectra cannot be used for spectral simu¬ 
lations at temperatures other than the one of the original 
experiment. Spectra l data on HNO 3 is o nly given in the 
HITRAN database llRothman et all l2013l l for up to 1770 
cm~^ (wavelengths longer than 5.6 pm). At shorter wave¬ 
length the Pacific Northwest National Laboratory (PNNL) 
database dSharpe et all l2004l i provides infrared cross sec¬ 
tions, but again these are only valid at the temperature for 
which they are recorded. Our results compare favourably 
with these sources, al though this comparison shows both 
have their limitations (IPavlvuchko et aLll2015all . 


The ExoMol project aims at generating comprehen¬ 
sive line observation and modelling atmospheres of ex¬ 
oplanets and other hot astronomical objects such as 
brown dwarfs and cool stars; its ai ms, scope and method¬ 
ology have been summarised by iTennvson fc Yurchenko! 
(l2012li . The project has provided rotation-vibration line 
lists f or a number of polyat omic molecules su ch as 
HCN (iBarber et al.ll2014|j . HaS (lAzzam et all l2015lL PH 3 
dSousa-Si lva et all 2015lL HaCO dAl-Refaie et al. 2015ll an d 


CH 4 ( Yurchenko fc Tennvsonll2014l : lYurchenko et al.llAlll j. 

None of these molecules contains more than two heavy atoms 
and the only pentatomic molecule considered, methane, has 
four hydrogens and high symmetry. It is clear that com¬ 
puting a comprehensive, temperature-dependent line list 
for HNO 3 represents a considerable computational chal¬ 
lenge. To address this challenge we have developed a hy¬ 
brid variational-perturbation theory procedure for com put- 
ing spectra of such molecules dPavlvuchko et aDl2015bl l and 
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have tested this for room -temperature nitric acid spectra 
dPavlvuchko et ^l2015al l. In this work we present a com¬ 
prehensive line list for nitric acid which should be valid for 
temperatures up to 500 K. Given that nitric acid is only 
likely to exist in atmospheres which also contain water, we 
explicitly consider regions where HNO 3 spectra are likely to 
be observable in a humid atmosphere. 


2 METHOD 



Rotation-vibration line li sts were generated using 
the program ANGMO L (iGribov &: Pavlvuchkol 1 19981 : 
iPavlvuchko et al]|2015lJ) . Although angmol is designed to 
solve the nuclear motion problem for polyatomic molecules 
using a hybrid variational-perturbation theory method, it 
actually provides a complete environment for performing 
the line list calculations. ANGMOL automatically generates 
the inputs required to drive the appropriate ab initio 
electronic structure progr am, here the qua ntum chem¬ 
istry package MOLPRO (I Werner et al.l [2012fl . to provide 
the necessary inputs for the nuclear motion calculations 
(potential energy, Hessian matrix and dipole moments). 
It also automatically adjusts both the potential energy 
surface (PES) and dipole moment function (DMF), which 
are represented as Taylor-expansions about equilibrium, 
to reproduce observed line pos itions and intensities using 
the method of regularisat ion llTikhonov fc Arsenin! 1 19791 : 
iGribov fc Dementiev! [l99^ ~l. This method uses constraints 
provided by the initial ab initio calculations to allow many 
more parameters to be varied than there are experimental 
data. For full det a ils of t he hts and associated surfaces see 
IPavlvuchko et ahl (l2015al l. 

ANGMOL solves a Watson-like nuclear motion Hamil¬ 
tonian expressed in internal curvilinear vibrational coordi¬ 
nates in three steps. First the rotation-less (J = 0) vibra¬ 
tional problem is solved in a basis of Morse oscillators for 
the stretches and harmonic oscillators for the bends by direct 
diagonalisation of a Hamiltonian matrix for low-lying vibra¬ 
tional states with contributions from higher states (basis 
functions) included using perturbation theory. In the sec¬ 
ond step, rotational problems for each vibrational state are 
solved using the eigenvectors from the vibrational problem 
and Wigner rotation matrices. An Eckart embedding is used 
to minimise Coriolis interactions which are again introduced 
using perturbation theory. The final step computes transi¬ 
tion intensities; this step reduces the size of the calculation 
by predicting which transitions will be too weak to make 
a signihcant contribution and not calculating them. Transi¬ 
tions whose intensity is less than 10 “^^ km/mole (« 10 “^^ 
cm/molecule.) were neglected. Integrals over the kinetic en¬ 
ergy operator, PES and DMF are all simplified by the use 
of Taylor expansions for the operators. Here the PES and 
DMF were represented as fourth-order and second-order ex¬ 
pansions, res pectively. Further details of the calculation can 
be found in (|Pavlvuchko et ahll2ni5ah . 

Our aim was to compute a HNO 3 line list which is com¬ 
plete up to 7000 cm“^ (longwards of 1.42 /rm) for tempera¬ 
tures up to 500 K. For this we considered energy levels up to 
9000 cm“^ and rotational states with J up to 100. J = 100 
corresponds to a rotational excitation of about 2130 cm~^. 

The initial vibrational (J = 0) problem was solved us- 


Wavelength {faro) 

Figure 1. Comparison of our compute d spectrum (Calc) with 
cross sections from the PNNL database dSharpe et al.ll2004^ at a 
temperature of 298 K. Cross sections were generated from our line 
list using a Voigt profile with a Doppler and pressure broadening 
widths both set to 0.075 cm~^. 


ing vibrational basis sets with polyad number 14. By this 
we mean that all combinations of functions were included 
when the sum of the orders of the polynomials was less than 
or equal to 14. The size of the basis required increases com- 
binatorially with polyad number N-. for HNO 3 , which has 9 
vibrational modes, the size of the basis is given by the bino¬ 
mial coefficient jv-i-gC'g which equals 817 190 for A" = 14 and 
48 620 for N — 9. Matrices constructed with polyad 9 were 
explicitly diagonalised once the effects of higher states were 
include d using perturbation theory, see IPavlvuchko et ahl 
ll2015ah for further details. These calculations gave a total 
of 22 049 vibrational states below 9000 cm“^. 

Rotationally excited states were computed based on cal¬ 
culation which only considered up to polyad 7 but with the 
vibrational band origins fixed by the larger J = 0 calcula¬ 
tion. This gave 9477 states ly ing below 7000 cm~^ of w hich 
intensity considerations (see IPavlvuchko et ahl (l2015al l for 
details) showed that only 1715 need to be considered ex¬ 
plicitly. No transition whose intensity was less than 10~^^ 
cm/molecule was considered. 

Our final line list links 17 494 715 vibration-rotation 
energy levels with 6 722 136 109 transitions. 

Temperature-dependent partition functions for HNO 3 
were calculated by explicitly summing all the calculated en¬ 
ergy levels. T able [J compares ou r results with those used 
by HITRAN (iFischer et al.l l2003h . The agreement is very 
good; our results are systematically slightly larger by less 
than 0.2% at low temperatures rising to 0.9 % at 500 K. 
This increase is probably due to our better treatment of 
anharmonic effects, particularly in the large-amplitude, low- 
frequency OH torsional mode. A file containing the partition 
function 1 K steps for temperatures up to 500 K is given in 
supplementary material. 


3 RESULTS 

The line list contain almost 7 billion transitions. For 
compactness and ease of use, it is divided into sep¬ 
arate energy level and transitio ns file. This is done 
using standard ExoMol format (iTennvson et al.l l2013ll 
which is based on a method o r iginal ly developed for 
the BT2 line list llBarber et al.1 1 200611 . Extracts from 
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Table 1. Partition function, Q(T), HNO 3 , as a function of tem¬ 
perature. 


T 

K 

Q(T) 

This work 

Q(T) 

Fischer et al. ('2003') 

60 

15032.00 

15010 

110 

37418.02 

37374 

160 

67172.12 

67105 
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107006.55 

106880 
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310 
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238250 

360 
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Figure 3. Absorption spectrum of a nitric acid (regular set of 
lines) and water vapour (four strong lines) mixture in the 1330 - 
1340 cm“^ region. The spectrum is for 300 K with a concentration 
of 2 % nitric acid. 


TRAN (|Rothman et al.ll2013l) are given bv IPavlvuchko et alJ 
ll2015ah . 

Given that HNO 3 is being considered as possible 
biomarker and that any biomarker will probably have to 
be observed in a humid atmosphere, we consider a num¬ 
ber of mixed water - nitric acid spectra. Figure [5] gives an 
overview of such spectra which have been constructed us¬ 
ing water from HITRAN with 2 % HNO 3 added for 300 K 
and 500 K. As is well-known, water absorptions dominate 
much of the infrared. However, HNO 3 features are clearly 
visible at about 7.5 /rm an d 11.2 ^m. These are precisely 
where HNO 3 is clearly seen in o ur own atmosphere from 
space telecka fc De Mazierel[l996h . 

Figures[3]and[4]give detailed comparisons for the two re¬ 
gions in question. The spectrum in the 7.5 fj,m region shows 
a clear, regular spectral signature of HNO 3 but would re¬ 
quire both a significant fraction of HNO 3 present and high 
resolution for a successful detection. When viewed in the 
Earth’s atmosphere this region also shows strong absorption 
lines due to CO 2 , although there are not enough of these to 
mask the HNO 3 . The 11.2 /rm region is actually much more 
promising. The hgure is drawn for only a trace concentra¬ 
tion of HNO 3 (0.01 % of water), much less than is present 
in the Earth’s atmosphere. Furthermore there is a distinct 


Wavelength (/xm) 

Figure 2. Temperature-dependent absorption spectra of a ni¬ 
tric acid (smooth curve) and water vapour (spiky curve) mixture. 
Spectra are for 2 % nitric acid at 300 K (upper) and 500 K (lower). 
The water spectra show a strong window around 4 pm. 


the start of the files are given in Tables [2] and [3] The 
full line list can be downloaded from the CDS, via 
ftp://cdsarc.u-strasbg.fr/pub/cats/J/MNRAS/xxx/yy 
or http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/MNRAS//xxx/yy 
The line lists and partition function, as well as the ab¬ 
sorpt ion spectrum given in cross section format (iHill et al.l 
I2OI3I) . can all be obtained from there as well as at 


HNO 3 band centred about 11.25 pm and a feature 11.38 
pm which should be visible (possibly blended together) at 
lower resolution. This would appear to be the most promis¬ 
ing region for an HNO 3 detection. We note that there is one 
further HNO 3 band clearly visible in Fig. [2] at about 3 pm. 
This feature is a blend of the 2 v 3 overtone (band centre at 
3404.4 cm~^), and the 1/2 + V 4 , (2998.5 cm“^) and 122 -f 1^3 
(3022.1 cm“^) combibation bands. This compound band sits 
in a window in the water spectrum but is not particularly 
prominent in the spectrum of the Earth’s atmosphere, pre¬ 
sumably because it is actually signihcantly weaker than the 
features at 11.2 pm and 7.5 pm. 


4 CONCLUSIONS 


WWW.exomol.com 

Figure [T] compares our calculation s with the cross sec - 
tions given in the PNNL database dSharpe et al.l l2004h . 
More detailed comparisons with both PNNL and HI- 


As part of the ExoMol project we have computed a com¬ 
prehensive line list for the possible biomarker nitric acid 
applicable to temperatures up to 500 K. To do this we 
have had to develop a novel methodology capable of treat- 
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Table 2. Extract from the state file for HNO 3 . Full tables are available from http: //cdsarc.u-strasbg.fr/cgi-bin/VizieR?-source=J/MNRAS/xxx/yy 


n 

E 

9 

J 

r 

VI 

V2 

VZ 

U4 

U5 

VQ 

V7 

VQ 

vg 

1 

0.000000 

6 

0 

A’ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

458.200000 

6 

0 

A” 

0 

0 

0 

0 

0 

0 

0 

0 

1 

3 

580.300000 

6 

0 

A’ 

0 

0 

0 

0 

0 

0 

0 

1 

0 

4 

646.450000 

6 

0 

A’ 

0 

0 

0 

0 

0 

0 

1 

0 

0 

5 

763.100000 

6 

0 

A” 

0 

0 

0 

0 

0 

1 

0 

0 

0 

6 

879.050000 

6 

0 

A’ 

0 

0 

0 

0 

1 

0 

0 

0 

0 

7 

896.300000 

6 

0 

A’ 

0 

0 

0 

0 

0 

0 

0 

0 

2 

8 

1038.000000 

6 

0 

A” 

0 

0 

0 

0 

0 

0 

0 

1 

1 

9 

1100.800000 

6 

0 

A” 

0 

0 

0 

0 

0 

0 

1 

0 

1 

10 

1148.975745 

6 

0 

A’ 

0 

0 

0 

0 

0 

0 

0 

2 

0 

11 

1205.600000 

6 

0 

A’ 

0 

0 

0 

0 

0 

1 

0 

0 

1 

12 

1213.857206 

6 

0 

A’ 

0 

0 

0 

0 

0 

0 

1 

1 

0 

13 

1275.720926 

6 

0 

A’ 

0 

0 

0 

0 

1 

0 

1 

0 

0 

14 

1289.000000 

6 

0 

A” 

0 

0 

0 

0 

0 

0 

0 

0 

3 

15 

1303.100000 

6 

0 

A’ 

0 

0 

0 

1 

0 

0 

0 

0 

0 

16 

1325.650000 

6 

0 

A’ 

0 

0 

1 

0 

0 

0 

0 

0 

0 

17 

1340.608559 

6 

0 

A” 

0 

0 

0 

0 

0 

1 

0 

1 

0 

18 

1343.600000 

6 

0 

A” 

0 

0 

0 

0 

1 

0 

0 

0 

1 

19 

1403.878302 

6 

0 

A” 

0 

0 

0 

0 

0 

1 

1 

0 

0 

20 

1450.282008 

6 

0 

A’ 

0 

0 

0 

0 

1 

0 

0 

1 

0 


n: State counting number. 

E\ State energy in cm~^. 
g: State degeneracy. 

J: Rotational quantum number. 

F: Total symmetry. 

ui, U 2 , U 3 , U 4 , U 5 , VQ, U 7 , ug, vg: vibrational quantum numbers. 



Wavelength (/2m) 


Figure 4. Absorption spectrum of a nitric acid (regular set of 
lines) and water vapour (isolated lines) mixture in the 875 - 900 
cm”^. The spectrum is for 300 K with a trace concentration of 
0.01 % nitric acid. 


ing the anharmon i c moti ons of a heavy 5-atom molecule 
(|Pavlvuchko et ^ ] l2015bll3l . The line lists, which con¬ 
tains almost 7 billion lines, can be downloaded from the 
CDS, via ftp://cdsarc.u-strasbg.fr/pub/cats/J/MNRAS/, 
or http://cdsarc.u-strasbg.fr/viz-bin/qcat7J/MNRAS/, or 
from www.exomol.com. 

Simulated spectra for a water-rich atmosphere suggest 
that there is real prospect of detecting trace quantities of 
nitric acid in the 11.25 /rm region. 


Table 3. Extracts from the transitions file 

for HNO3. Full tables are available from 
http: //cdsarc.u-strasbg.fr/cgi-bin/VizieR?-source=J/MNRAS/xxx/yy 


F 

I 

Api 

14516084 

14516083 

4.5130E-10. 

14516265 

14516264 

8.1767E-10. 

14515899 

14515898 

1.8398E-09. 

14516082 

14516081 

2.4872E-08. 

14516263 

14516262 

2.8596E-08. 

14515897 

14515896 

7.7843E-08. 

14516080 

14516079 

6.5506E-07. 

14516261 

14516260 

6.2021E-07. 

14515895 

14515894 

1.5810E-06. 

14516078 

14516077 

7.7309E-06. 

14516080 

14516078 

7.1042E-06. 

14515895 

14515893 

1.5567E-05. 

14516259 

14516258 

5.9759E-06. 

14515897 

14515895 

1.8498E-05. 

14515893 

14515892 

1.4427E-05. 

14516261 

14516259 

1.2759E-05. 


F: Upper state counting number; 
I: Lower state counting number; 
Apj: Einstein-A coefficient in s“^. 
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